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Abstract 

Whenever real particle production occurs in quantum field theory, the imag- 
inary part of the Hadamard Elementary function is non- vanishing. A 
method is presented whereby the imaginary part of may be calculated 
for a charged scalar field in a static spherically symmetric spacetime with 
arbitrary curvature coupling and a classical electromagnetic field . The 
calculations are performed in Euclidean space where the Hadamard Elemen- 
tary function and the Euclidean Green function are related by \G^ = Ge- 
This method uses a 4 th order WKB approximation for the Euclideanized mode 
functions for the quantum field. The mode sums and integrals that appear in 
the vacuum expectation values may be evaluated analytically by taking the 
large mass limit of the quantum field. This results in an asymptotic expansion 
for in inverse powers of the mass m of the quantum field. Renormaliza- 
tion is achieved by subtracting off the terms in the expansion proportional to 
nonnegative powers of m , leaving a finite remainder known as the "DeWitt- 
Schwinger approximation." The DeWitt-Schwinger approximation for G^ 
presented here has terms proportional to both m~ l and m~ 2 . The term pro- 
portional to m~ 2 will be shown to be identical to the expression obtained 
from the m~ 2 term in the generalized DeWitt-Schwinger point-splitting ex- 
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pansion for The new information obtained with this present method is 

the DeWitt-Schwinger approximation for the imaginary part of , which is 
proportional to m _1 in the DeWitt-Schwinger approximation for derived 
in this paper. 
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I. INTRODUCTION 



Despite the abscence of a full quantum theory of gravity, there are well known methods 
to calculate useful quantities involving quantum fields in curved space. One of the most 
successful of these methods is the semiclassical approach, wherein the gravitational and 
electromagnetic fields are treated classically while the various matter fields that act as sources 
are treated quantum mechanically. In this approach, the evolution of the spacetime from a 
set of given initial conditions is described by the combined Einstein-Maxwell Field Equations, 

G iu , = 8tt(T iu/ ), (1) 

and 

F^ = 4tt(^), (2) 

where (T^ u ) and are the vacuum expectation values (VEVs) of the stress-energy ten- 
sor and the current due to a charged quantized scalar field. This paper follows the sign 
conventions of Misner, Thorne, and Wheeler |L] and uses natural units {G = c = % = 1) 
throughout. 

In the semiclassical regime, the calculation of quantities such as (</> 2 ) and (T^) may 
begin with calculation of the Feynman Green function Gp(x,x') defined by 

tG F (x, x) =< O(T[0(x)0(x')])O >, (3) 

where T[. . .] is the time ordering operator, <p{x) is the quantum field evaluated at the sta- 
tionary spacetime point x, and 4>(x') is the quantum field evaluated at the nearby point x'. 
One of the major difficulties in using the semiclassical method in quantum field calculations 
is the fact that infinities appear in the VEVs in Eqs.([l])-(||) in both flat and curved 
|2||| spacetimes. Various methods have been used to isolate those infinities and remove them 
from the physical theory ||, yet arguably the most powerful of these methods is based on 
the classic work of Schwinger 0]. 
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Schwinger calculated the Feynman Green function for a quantized fermion field by intro- 
ducing a fictitious, non-quantum mechanical Hilbert space. This (4 + l)-dimensional Hilbert 
space was constructed from the 4 spacetime dimensions in addition to a fifth dimension which 
was identified as the proper time parameter s within this fictitious space. Working within 
this proper time space, Schwinger was able to both isolate the divergences in the quantum 
field integrals for the fermion current (j^{x)) = \im x i^ x ietr^Gp^x, x')] produced by an 
external electromagnetic field in flat spacetime, and to use those divergences to renormalize 
the charge e of the quantum field. 

Schwinger's calculations of the quantum action functional W were performed by trans- 
forming the integrals involving the 4 spacetime dimensions into the momentum represen- 
tation. The integrals were evaluated using perturbation expansions in powers of eA^ and 
eF^ u , where and F^ u are the gauge field vector and the electromagnetic field tensor, 
respectively. Were these integrals to be evaluated without further modification, Schwinger 
noted H that conservation of energy and momentum would dictate that no pair creation 
would occur, or = 0. Schwinger's solution to this situation was to add an infinitesimal 
imaginary part to the action integral, resulting in W acquiring a positive imaginary part. 
Schwinger associated this imaginary part of W with pair production by stating that 
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e -2Im[W] ^ 



represented the probability that no pair creation would occur. After this identification was 
made, Schwinger's calculation of the pair creation rate yielded a series expansion for the 
probability of pair creation per unit four volume T by a constant external electric field, 

r = ^Y.n- 2 e^), (5) 

^" n=l 

where m and e are the mass and charge of the fermion field, respectively, and E is the 
magnitude of the electric field. The expansion of Eq.(|3p yields a series involving inverse 
powers of m 2 . 

Schwinger's calculation for the pair creation rate assumed that spacetime was flat. Build- 
ing on the work of Schwinger, DeWitt |||| extended the proper time method to include 



curved spacetime results. Using what is now called the DeWitt-Schwinger proper time 
method, DeWitt calculated an asymptotic expansion for the Feynman Green function for a 
real scalar field in an arbitrary curved spacetime. This asymptotic expansion was an expan- 
sion in inverse powers of m 2 , where m is the mass of the quantum field. This expansion of 
Gp(x, x') yielded an expression with both real and imaginary parts. The real and imaginary 
parts of Gf(x,x') could then be identified with other Green functions by using the identity 

Gf(x,x') = G(x,x') — -iG^\x,x'), (6) 

where G(x, x') is one-half the sum of the advanced and retarded Green functions, and the 
Hadamard Elementary function G^(x,x') @ is defined by 

GW(x,x') = (O|{0(x),0(x')}|O). (7) 

Thus, the DeWitt-Schwinger asymptotic series expansion for G^'(x, x') has the general form, 

G {1) {x,x r ) ~ B +2 (x,x')m 2 + B {x,x')m° + B„ 2 (x,x')m- 2 + . . . , (8) 

where the B n are coefficients constructed from curvature tensors. 

In the semiclassical method, quantities such as ((f) 2 ) and (T^) may be calculated from 
G^ (x, x') by taking the appropriate derivatives with respect to either the stationary space- 
time point x or the nearby point x' |6|,p|-p|,31 ; e.g., 

(<?) = hmV) (*,*'), (9) 

x'—>x Z 

and 

(T^(x)) = ]m jfe[i(±-£)(^G^+^^ (10) 

In Eq.(|T0"D, the vertical bar "|" indicates gauge covariant differentiation, and g'V is the 
bivector of parallel transport ||, which is used to transport the information from the nearby 
point x' to the stationary point x. The infinities appearing in the VEV in Eq.(^) lead 
to the well known divergences in the unrenormalized expression for (T^^x)) in Eq.(|I"0"D 



@. DeWitt's asymptotic expansion of G^(x,x') in powers of m 2 isolated the infinities, 
or the "counterterms," of both G^(x,x') and (T^), with the infinities being found as 
coefficients of the m 2 and m° terms in the series of Eq.(H) |T0f] . Once these infinities have 
been subtracted from the unrenormalized VEVs in Eqs.(p|)-([T0|), then the finite remainder 



carries the information about the physics of the spacetime. While these subtractions may 
be performed in principle, the remaining expressions are often so complicated that they can 
not be evaluated analytically 0. 

Investigators have used the semiclassical method to calculate renormalized expressions for 



both the vacuum polarization (<jr) ren |l2|-^7[, and the stress energy tensor (T^ u ) ren [|T8| -|2T 



in specific spacetimes. In each of these cases, the point-splitting counterterms were used to 
renormalize the divergent VEVs. These calculations were performed in spacetimes whose 
high degrees of symmetry allowed the subtractions of the counterterms to be evaluated 
analytically. See Ref. |2l[ for a more thorough discussion of the difficulties encountered in 
these subtractions. 

The DeWitt-Schwinger point-splitting expansion for may also do more than simply 
isolating the divergent counterterms necessary to renormalize the VEVs encountered in 
semiclassical quantum field theory. If the terms proportional to nonnegative powers of m are 
discarded, the remaining terms in Eq. (§) constitute the "DeWitt-Schwinger approximation" 
for for fields whose masses are large when compared to the magnitude of the 



coefficients B n . Retaining the first one or two of these terms has been shown to provide a 
close approximation to the exact results in several cases P~^JTP|jr7|f2T[ ] . 
In the case of charged scalar fields, G^ involves the complex quantity 

G«(x,x') = (O|{0(x),0*(x')}|O), (11) 

where (j>(x) is the charged scalar field. Even though G^ is generally complex, since all of the 
coefficients B n are real, Eq.(|8|) only contains part of the information about G^. DeWitt has 
pointed out that the point-splitting asymptotic expansion of G^ in inverse powers of m 2 
is incapable of yielding the imaginary part of G^\ Yet whenever real particle production 



occurs (jW will have a non-zero imaginary part. This is readily seen by examining the 
expression for the current due to a charged scalar field, 

I IP r i 

f = = 2 [ {jD ^' ^ ~ {D "^ ' (12) 

where e is the charge of the field, = (V M — ieA^) is the gauge covariant derivative, 
is the background gauge field, and the asterisk * denotes the complex conjugate. Using 
Eq. (|TTD , and making the transition from classical to quantum fields ||24|| , gives the vacuum 



expectation value of the current due to the coupling between the charged scalar field and 
the background gauge field, 

(f(x)) = lim % \(G^ + g^G^') - (g^ + g^G^'Y] . (13) 

x'^x 4 ^ V / \ / - 

Since the coefficients B n in Eq.(§) are real, using the DeWitt-Schwinger point-splitting 
expansion involving only inverse powers of m 2 to construct the current (j^) from Eq.(|l3"D 
can only yield a value of zero. Yet the current is not identically zero, as evidenced by 
Schwinger's calculation in flat space, thus illustrating the limitation the DeWitt-Schwinger 
point-splitting expansion suffers by yielding only the imaginary part of G^ 1 ' . 

In this paper, a method is presented to obtain an analytic expression for the imaginary 
part of in a static spherically symmetric spacetime with a gauge field M 1 . The imaginary 
part arises in a straightforward manner due to the presence of the gauge field N 1 and does 
not require the addition of any extra terms during the course of the calculation. The 
renormalized expression GfJ n would ordinarily be constructed by subtracting the DeWitt- 
Schwinger point-splitting counterterms from the unrenormalized expression for G^; 

_ (7(1) _ (iA) 



The mode sums and integrals in Eq.(fbl|) can not be evaluated analytically at present due to 
the complexity of the terms in the infinite sums and integrals. A method recently described 
by Anderson, Hiscock, and Samuel (AHS) [21J yields a way to construct an approximate ex- 
pression for G$ n . All subtractions are performed in Euclidean space where the relationships 
between various Green functions are given by |~7| 



((f) 2 ) = ^G w {x, x') = tG F (x, x') = G E {x, x'). (15) 
Using the Euclideanized metric, the mode functions in the VEV of Eq . (|TTD may be rewritten 



in a WKB form JT^J2T|. The renormalized expression GE,ren is constructed by moving into 



Euclidean space and performing the subtraction 

GE,ren = Ge,WKB ~ Ge,WKB(Hv, (16) 

where Ge,wkb is the expression obtained in Euclidean space upon substituting the Euclidean 
WKB mode functions into Eq. ([TT|) , and Ge,wkb<Hv contains the ultraviolet divergences found 
in Ge,wkb- The mode sums and integrals in Eq. (|T6|) still can not be evaluated analytically 
in their present form. In order to construct an expression for G_E )ren , the fourth order 
WKB approximations for the exact mode functions are substituted into the quantum field 
expressions. Then, the mode sums and integrals in the subtraction of Eq.(|T6D may be 
evaluated analytically by expanding the integrands and summands in powers of the mass 
m of the quantum field. This mass is assumed to be large when compared to the inverse 
of the radius of curvature of the spacetime. By working within this large mass limit, a 
DeWitt-Schwinger asymptotic expansion results for Ge^u similar to Eq.(H); 

GE,ren = (Ge,WKB ~ G E,WKBdiv) large m (17) 

~ B +2 {x, x')m 2 + B + i(x, x')m x + B (x, x')m° (18) 
+B^ 1 (x 1 x')m~ 1 + B. 2 (x,x')m~ 2 + . . . . (19) 

Eq.(|19"D differs from Eq.(|8D in that the presence of the gauge field in this method is directly 
responsible for the presence of the new terms involving odd powers of m. Renormalization 
is achieved by discarding terms proportional to nonnegative powers of m, leaving a finite 
DeWitt-Schwinger approximation for GE,ren] 

G E .ren ~ B^(x, x')m~ l + B_ 2 (x, x')m~ 2 . (20) 

In Euclidean space, all of these terms are real. When the final expression is rotated back 
to the Lorentzian sector, the real and imaginary parts of G^ are obtained. The terms 
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proportional to m" 1 are all imaginary and are proportional to both odd powers of the gauge 
field A 11 and its derivatives along with odd powers of the charge e of the quantum field. The 
terms proportional to m~ 2 are all real and are exactly those obtained from the m~ 2 term of 
the DeWitt-Schwinger point-splitting expansion p!0|Ji7| , |22[] . 

In Sec. |IT], an expression is derived for the unrenormalized Euclidean Green function 
Ge in a static spherically symmetric spacetime. In this paper, the field is assumed to 
be at zero temperature, although this method would extend to the analysis of fields at 
nonzero temperature | 17jl . Sec. |i| renormalizes the expression for Ge derived in Sec. |J. 
The renormalization subtractions are evaluated analytically by evaluating the mode sums 
and integrals in the unrenormalized expression in the limit that the mass of the quantum 
field is large when compared to the inverse of the radius of curvature of the spacetime. 
The resulting DeWitt-Schwinger approximation will be shown to contain both the real and 
imaginary parts of G^ 1 ' , with the main goal of this paper being obtaining an expression for 
the imaginary part. 

Sec. [TV] contains the results of the calculation of the asymptotic expansion of G^\ The 
terms proportional to m~ 2 will be shown to give the same results as the m~ 2 term in 
the generalized DeWitt-Schwinger point-splitting expansion when the general expression is 
evaluated in a static spherically symmetric spacetime |10|JT7[| . The terms proportional to 
m _1 are all real in Euclidean space but are all imaginary when rotated back to Lorentzian 
space. All of the m~ l terms are proportional to odd powers of the charge e of the complex 
scalar field, as well as odd powers of the gauge field A 11 and derivatives, and thus these terms 
all vanish when either the quantum field is uncharged or the spacetime is uncharged. This 
connection to previous work will serve as a check on the validity of the method presented 
here. 
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II. UNRENORMALIZED GREEN FUNCTIONS FOR THE CHARGED SCALAR 

FIELD IN EUCLIDEAN SPACE 



The goal of this paper is to derive an expression for the imaginary part of so 
that curved space phenomena involving real particle production due to the presence of an 
electromagnetic field may be investigated. Thus it is useful to begin with the electromagnetic 
vector potential for a static spherically symmetric electric field like that of the Reissner- 
Nordstrom spacetime; 

A" = (A*(r), 0,0,0). (21) 



The calculations will be simplified by moving to Euclidean space ||17|| . The Euclidean metric 
for a static spherically symmetric spacetime is 

ds 2 = f(r)dr 2 + h{r)dr 2 + r 2 dtt 2 , (22) 

where r = —it, f(r) and h(r) are the same functions as in the Lorentzian metric for a static 
spherically symmetric spacetime 

ds 2 = -f{r)dt 2 + h{r)dr 2 + r 2 dQ 2 , (23) 

and t is the coordinate whose derivative (d/dt) is the timelike Killing vector in the Lorentzian 
sector [17|. The Lorentzian and Euclidean expressions for the invariant A^A^ are given by 



A^A* = guA* L A? L = -f{r)A* L A* L = f(r)(iAl)(iAi) 

= g TT AlA T E = f(r)A^ E , (24) 

where the subscripts "L" and "E" refer to the Lorentzian and Euclidean sectors, respectively. 
Eq.(p4|) indicates that 



A t e = Ml- (25) 
This shows the relation between the gauge invariant derivative operators in both sectors is 

10 



g^DlDl = g tt (V - zeA*)(V* - ieA 1 ) + ■■■ 

= -/(r)(V* - 2eA*)(V* - zeA*) + • • • 
= /(r)(iV* - ze(zA*))(zV* - ie(iA*)) + • • • 
= /(^)(V r - zeA T )(V r - ieA T ) + • • • 
= g TT D T D T + ■■■ 

= g, v D E D E (26) 

indicating 

£> T = V r - ieA T (27) 

for the Euclidean sector. 

An explicit form for the Euclidean Green function for fields at zero temperature is needed. 
The case of zero temperature is chosen here for simplicity, and the derivation here follows 
those of Anderson [I7J and AHS [^TJ. The action of the wave operator on the Euclidean 
Green function is given by 

[g, u D E D E - (m 2 + £R)]G E (x, x') = -g- 1/2 S\x, x') 

5(T,T')5(r,r>),5(n,n>) 

r 2(f h y/2 ■ y ) 

The delta function 5(Q, Q') may be expanded using the Legendre polynomials with the result 

1 oo 

S(n, 0) = — ^(2Z + l)P,[cos( 7 )], (29) 

where cos(7) = cos(6 l ) cos(6 |/ ) + sm(9) sm(9') cos(0 — </>'). With the scalar field chosen to be 
at zero temperature, the delta function for the split in the r — r' direction is given by 

$(t,t') = ±- r duo e^'\ (30) 

Z7T J-oo 

Thus, the Euclidean Green function for a field at zero temperature is 

1 /•oo 00 

G E (x,x') = —J o duo cos^(r-r')]E(^+l)^[cos(7)]^(r,r'), (31) 



1=0 
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where S w i(r,r') is an unknown function of the coordinate r. 

Applying the wave operator of Eq.(|26|) to Gg gives a differential equation for S^r, r'), 



ld 2 S 



h dr 2 



+ 



2 1 df 1 dh] dS 
rh 2fh dr 2h 2 dr dr 

(u-^iv + v^ R 



f 



s 



S(r, r') 



(32) 



where m is the mass of the field, and A T is the only non-zero component of the Euclideanized 
vector potential. S^r, r') is given by 



S u i(r,r') = C u ,ip ul (r)q u i(r f ), 



(33) 



where C w i is a normalization factor, and Pu,i{r) and q w i{r') are the solutions of the homoge- 



neous form of Eq. (|32|) . In Refs. [17j and [ 2lj , S^^r, r') is written in the form 



Sui(r,r') = C w ^(r < )g^(r > ). 



(34) 



S^iir, r') satisfies the Wronskian condition 



a 



Ull 



Pur 



dqcoi dp u 

-j Qui —r~ 

dr dr 



1 (K 



1/2 



(35) 



obtained by integrating Eq.(|32|) once with respect to r from r — e to r + e, with e — > in 
the end. 



The functions p w \ and q w i may be put into a WKB form using ||17 

1 



J r W(r')(!f) 1/2 dr' 



Puji{r) 



2r 2 W(r) 

1 p -f W (r>)($y /2 dr> 



(36) 

(37) 



l2r 2 W{r) 

In Eq. fl36|) , Puji(r) is well behaved at r = and the event horizons, but is divergent at 
r = oo. The function qu( r ) * s divergent at r = and the event horizons, but is well behaved 
at r = oo. See Anderson |T7J for a more complete discussion of the properties of these 
functions. Eqs. (P5])-(P^) show that C w \ = 1. Using Eqs.(|3"ED and (|3~TD, Eq.(|3"2]) assumes the 
familiar WKB form 
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w 1 



n 2 (r) + V 1 (r) + V 2 (r) + 



f d 2 W (Idf 1 dh\ 1 dW 3/ / 1 dW' 



hW dr 2 \hdr fh 2 dr 2W dr 2 h \W dr 



(38) 



where 



Q 2 {r) 
V x {r) 
V 2 (r) 



= [u 



eA T f + m 2 f + (/ + 



1x2 / 



Id/ / d/t 



2r/t dr 2r/i 2 dr 

tRf = -e/ x 

1 d 2 / 1 /d/ 



4 r 2' 



(39) 
(40) 



1 dfdh 2 df 



+ 



2 d/i 2 
+ 



fhdr 2 2f 2 h \dr I 2fh 2 drdr rfhdr rh 2 dr r 2 h 



(41) 



and where the product /(/ + 1) has been factored into two pieces according to 

l(l + l) = (l + \) 2 -\. (42) 

The simple form of the vector potential in Eq. (|2l| ) causes the WKB equation to assume 
a new but straightforward form. Previous work had WKB equations just like Eqs.(|38"|)-(f4T|) 



with A T = (see the references in Ref. |2T| for a more detailed list of work by others). 
The new feature in the WKB differential equations is the presence of the electromagnetic 
vector potential. Eq.(|38"D indicates how the complex quantities necessary for determining the 
charged scalar field current (j^) arise in the present method. The WKB function W(r), while 
real in the Euclidean sector, is now complex in Lorentzian space due to the electromagnetic 
field of Eq.fl3|). 

Eq.Q3"8]) is not in general exactly solvable. Fortunately, it may be solved approximately 
by iteration. For example, the zeroth order solution for W is W = Q. Substituting W — > Q 
into the right hand side of Eq. d38|) and solving for W yields the second order WKB solution 

w< 2 >. 



W V = S l + ±(V 1 + V 2 ) + ±.V? 



1 

+ 4 



f d 2 n fidf _ f dh\ i dn _ 3 / i /d^ s 

W 2 !^ 2 + \hdr~ ~ h 2 ^) 2^2 dr" ~ 2hQ? I dr 



(43) 
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The fourth order solution, W^, is obtained by making the substitution W — > in 
the right hand side of Eq.(^8|) and solving for W. Just as in point-splitting, this iteration 
procedure could continue indefinitely. However, in the present case of obtaining the -^(m -1 ) 
and ?9(m -2 ) terms in the DeWitt-Schwinger expansion for the iterations need only be 
continued until the fourth order WKB solution, is reached. As pointed out by AHS 

2TJ] , a second order WKB expansion contains the information of the d(m 2 ) and d(m°) terms 



in the DeWitt-Schwinger point-splitting expansion, while a fourth order WKB expansion 
contains the information of the d{m~ 2 ) terms. Thus, the infinite counterterms in the theory 
may be reproduced by performing the present calculation using only W^ 2 \ while the DeWitt- 
Schwinger approximations of the t?(m -1 ) and $(m -2 ) terms require the use of W^ A \ 

For purposes of simplicity and to make connections to previous work |21|| , the choice is 
made to split the points along the r — r' direction such that (x, x') — > (r, r, 9, <fi; r', r, 9, (p). 
Making the definition 

A 1= lim£ 2{l + -) Pull {r)qM) =E , ( 44 ) 

r ^ r l=t L 2 J ^ L r 2 W _ 

with the limit r' — > r being taken in the last expression, allows the unrenormalized Euclidean 
Green function to be written as 

1 r°° 

G E ,unren(x,T;X,T') = {(j)) U nren = 77 / duj COs[u(t - t'))Ai, (45) 

47T JO 

for the separation along r — r'. The next section discusses the identification and removal of 
the two types of divergences that appear in this unrenormalized expression for Ge{x, t; x, t'). 



III. RENORMALIZATION OF THE A™ ORDER WKB APPROXIMATION IN 

THE LARGE MASS LIMIT 

There are two types of divergences which appear in Eq.(|45|). The first of these is a 
superficial divergence which appears when evaluating the sum over / with an upper limit of 
I = 00, while the second is an ultraviolet divergence due to performing the integrations over 
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u) with an upper limit of u = oo. This section presents the details of how the subtractions of 
these divergent terms from the unrenormalized Green function may be performed and how 
the resulting expression may be evaluated analytically. 

As discussed more fully elsewhere pl| . |i7|JT5|JT8|] , the superficial divergence due to the 
sums over I can not be real since Ge, unren must remain finite with the points separated as 
they are in Eq.(^). The sum and integral present in Eq. (f45|) will be evaluated later by 
assuming the mass m of the quantum field is large enough to expand the summand of A\ 
in inverse powers of m. At present, it is useful to consider how the false divergence over / 
is isolated assuming the sum and integral of Eq . (f45|) were to be evaluated without such a 
large mass expansion. 

The superficial divergence that appears in A\ may be identified by evaluating the sum 

OO 1 OO (111) 

£2(Z + i)^ = £^ (46) 

1=0 Z 1=0 ' vv 

in the limit I — > oo. The only part of W that contributes as I — ► oo is the zeroth order term 
of W, or Q; terms in W proportional to negative powers of Q do not diverge as I — > oo. 
Substituting W = fl into (/ + \)/r 2 W and expanding in inverse powers of Z, and then 
truncating the expansion at gives 

~ 2(1+1) ^ 1 ,AT\ 

> hm - 2 = > — r . (47) 

Thus, the counterterm l/(r/^) must be subtracted from the summand of Eq.([B|) in order to 
remove the superficial divergence over /. This now indicates the unrenormalized expression 
for Ge(x, r; x, r') is given by 

■(/ + §) 1 



r 2 W r/3 



(48) 



G E ,unren(x, T\ X, t') = / du COs[u;(r - /)] £ 

47T JO l=Q 

This technique will serve to remove any superficial divergences over I that may appear. In 
Eq. (^8|) , only the term l/(r/a) is subtracted from the single summand (I + ^)/r 2 W. Later, 
when the WKB approximation for W is substituted into Eq. (fS>D and the summands and 
integrands are expanded in the large mass limit, many terms involving factors of / will arise. 
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Only some of those terms will contain a divergence in /, and those that do must have the 
appropriate counterterms in I subtracted from them to remove this superficial divergence. 

Eq.pSD also contains the ultraviolet divergences that are known to arise in semiclassical 
quantum field theory. These divergences occur as u — > oo and must be subtracted from the 
unrenormalized expression. As in the case of the divergence over I, these subtractions must 
be performed term by term after the large mass expansions of A\ have been performed. 
Thus, the generalized form for the renormalized expression for Ge is given by 
1 f°° 

G E ,ren = ~. q / du COs[u(t - T )] X 
47T JO 

(49) 



47T 2 JO 

< — — - > — (/ — > 00 terms) — (u — > 00 terms) 



. r 2 W , 

=0 / WKB expansion terms 

where the subscript "WKB expansion terms" indicates all of the terms that result from a 
4 th order WKB expansion of (I + ^)/r 2 W, and the terms identified as divergent in I and u 
are subtracted from each of these terms of the expansion. 

The calculation of all the terms in Eq.(pi) is a straightforward but tedious process. The 



general procedure follows Appendices D, F, and G in the recent work of AHS |2T[ and will 
be described here. The presence of the gauge field increases the difficulty of the calculations 
substantially but the method of AHS can be modified to account for the gauge field. 

The calculation of G_B,ren starts with the substitution of the 4 th order WKB approxima- 
tion for W into the definition of Ax and the subsequent identification of the various I — > 00 
counterterms. To keep track of the order in the WKB expansions of quantities it is useful to 
introduce the dimensionless parameter a, letting a — > 1 in the end. In Eq. (p8|) , Q is i?(a°), 
Vi(r) is ^(a 1 ), and the rest of the terms are ?9(a 2 ). The quantity (I + \)/{r 2 W) is expanded 
in powers of the WKB parameter a, truncating the expansions at $(a 4 ), and setting a = 1. 
This will result in over 600 terms whose general form is given by 



J hjk — 

1=0 



2 ■ 



— (/ — ► 00 terms) 



(50) 



where f2(r) is defined in Eq.(|59"D, and various factors of f(r), Vx(r), eA T (r), . . . , that are 
present in front of each of the Lhjk have been omitted. This definition of Lhjk is similar to 
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that defined as Ljk in Ref. [pljj . In the present work, the presence of the gauge field has led 
to the need to keep track of additional factors of oo which now enter into the definition of 
the Lhjk- The notation here is chosen so that it most closely matches the notation of AHS. 
Determining the I — > oo subtraction terms proceeds as before by expanding the summands 
in each Lhjk in invervse powers of I and truncating the expansion at $(Z°). For example, the 
expression for L oi requires one subtraction term; 



>ooi 



E 

1=0 



2(1 + 



n 



2r 



(51) 



Determining the oo — > oo subtraction terms requires using the Plana sum formula |2T| . |25 
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[g(k + it) - g(k - it)] . 



(52) 



The first two terms of Eq.(|52"D may be computed analytically, while the third term is not in 
general able to be computed analytically. Since the ultraviolet, or u — >• oo, behavior of these 
sums is dominant, the third term of Eq.(|52|) may be expanded in inverse powers of u. The 
expansions are truncated at d(uj~ l ) since terms of $(u~ 3 ) are not ultraviolet divergent pT| . 
Then, the integrals that arise from the expansion in inverse powers of u may be computed 
analytically. Applying this procedure to L 00 i, for example, gives 
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(53) 



/ / V12 

The only for which there will be oo subtraction terms are those for which {j, k} have the 
values {0, 1}, {0, 3}, {1, 3}, {1, 5}, {2, 5} and {2, 7}; only h = occurs for these combinations 
of {j, k}. Comparison of L oi in this example with L i of AHS, along with Eq . (|25"D , show how 
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the gauge field will bring complex quantities into the expressions for the Lhjk and eventually 
into 

Once the terms divergent in both / and u have been determined for each of the Lhjk, the 
expression for GE,ren is now given by 

1 f°° 

G E ,ren = ^ J Q duJ C S K r ~ T ')]{ F ( r )( L hjk ~ ^/yfc,w-»oo)}, (54) 

where each of the F(r) above represent the functions such as f(r), Vi(r), eA T (r), . . . , that 
appear in the calculation but are not affected by the summations and integrations, and the 
curly brackets {. . .} indicate the large number of terms arising due to the WKB expansions. 
In Eq. (|54j) , the uj — > oo counterterms for each L h jk are indicated symbolically by Lhjk,u-*oo- 
Since the terms which diverge as u — > oo are subtracted in Eq.([54]), then the function 
cos[u;(t — t')] will be dominated by the infinitesimal separation (r — r'). This indicates 
that cos[u;(r — r')] ps 1 is valid. For the purposes of completeness, it should be noted that 
future work involving both and (T^) will involve more factors of u than those that 
arise within the Lhjk- Thus, it is useful to make the definition 

1 f°° 

Sihjk = -^—^ du La 1 (Lhjk ~ Lhjk,w-^oo) i (55) 

where the indices ihjk have been chosen so the present work is more easily compared to that 
of AHS. In the calculation of the Green function in this paper, % = throughout; in AHS, 
i = for the Green function calculation while i = or 2 when calculating (T^). The sums 
and integrals in Eq.fl55D still can not be evaluated analytically in general in their present 
form. They may be put into a form that can be evaluated analytically by taking the large 
mass limit of the quantum field. 

Each of the S^k may be computed in the large mass limit by first expanding each 
vfl (I + |) 1+2j /fi fc term in the Lhjk using the Plana sum formula of Eq. (|52|) . The first 
two terms of the Plana formula may be computed exactly. Again, the third term may 
not be computed analytically, but once expanded in the large mass limit each piece of 
the expansion may be computed analytically. The large mass limit of this third term is 
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obtained by expanding the integrand of the third term in inverse powers of the large quantity 



[uj — eA T ) 2 + m 2 f + (I + \) 2 f/r 2 . Each integral may then be evaluated analytically. After 
the integrations, the entire expression for each Sihjk is expanded in inverse powers of m. In 



Appendix G of Ref. ||21|| , AHS give an example of this type of calculation when they calculate 
the large mass expansion of Sooi- Here, the large mass expansion of Soooi w ^ °e shown for 
comparison. 

The / — > oo limit for Looi is found in Eq . (|5l|) , while the u — > oo subtraction term is 
found in Eq. (|53|) . This gives the expression to be evaluated in the large mass limit as 
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Using the Plana sum formula and expanding in inverse powers of the large quantity 



'uj - eA T ) 2 + m 2 f + (l + \) 2 f/r 2 yields 
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+ 105r 4 [(w - eA T ) 2 + m 2 f + ^fl 2 + " " " 

( 2r2 f 1 2 2\ 1 2er 2 A T \ ) 

cu+ — -mV - + ——- ) >. (58) 
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Performing the integrals over u and expanding in the large mass limit yields 
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64512m% 2 r 4 v ; v 7 

In Eq . (|59|) , the infrared cutoff parameter A has been introduced due to the lower limit u = 
on the integrals over u @J21|. Comparison of Eq.© with Eq.(G2) of Ref. |l| shows that 



the presence of the gauge field has brought odd powers of m into the asymptotic expansions 
for the Sihjk- Fortunately, when the large mass expansions of all of the S^jk are computed, 
many are found to be of such a high order in inverse powers of m that they do not contribute 
to the final expression for GE,ren since the DeWitt-Schwinger expansions are truncated at 
$(m~ 2 ). Note that the gauge field has also contributed terms of ^(m 1 ) and -^(m ) to Soooi, 
as it does for some of the other Sihjk- In this work, only terms which contain negative 
powers of m will be considered in the final DeWitt-Schwinger approximation. Finally, all of 
the S^jk are substituted into the general expression for GE,ren 

G E ,ren = {F(r)S ih jk}, (60) 

where the curly brackets {. . .} indicate the great number of terms to be added together, and 
like powers of m are collected. 
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IV. RESULTS 



When all the substitutions are made and terms of like powers of m are collected, the 
result is that the coefficients of the negative powers of m constitute the DeWitt-Schwinger 
approximation to Gs,ren for a complex scalar field in a static spherically symmetric space- 
time. A complex scalar field was chosen in order to make a connection to the earlier work of 
DeWitt and Christensen [O involving a real scalar field since the results of the present 
work must reduce to their results when the charge of the field vanishes. In their work, the 
generalized DeWitt-Schwinger approximation for the Hadamard elementary function for a 
scalar field is 
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where [a 2 ] is given by P,|2l 

H = 



— R al3 R af 3 + —R^Ra^s + -(- - R ;a 



2 V 6 ^ J 



12" 



:F af3 F aP . 



(62) 



In the present work, the m 2 term in the DeWitt-Schwinger approximation for Ge,; 
L > 2 ) is given in Euclidean space by 
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960 / h 3 480 fh 2 ) 

where the primes indicate differentiation with respect to r. This is exactly the same result 
that is obtained by writing Eqs.(|B"T|)-(|6"2"|) in terms of the metric functions of Eq.(^). Since 
the gauge field only appears here in the form (eA' T ) 2 , rotating back to Lorentzian space to 
obtain G^- 1 ' simply results in the appearance of a negative sign for this term. This demon- 
strates the correspondence between the present method in a static spherically symmetric 
spacetime and the method of the generalized DeWitt-Schwinger expansion; they yield the 
same information for the real part of in the limit that the mass of the quantum field is 
large when compared to the inverse of the radius of curvature of the spacetime. 

With the information in Eq . fl63|) previously known and accepted, it is fortunate that the 
present method not only duplicates that information but also yields more information. The 
DeWitt-Schwinger point-splitting expansion is fully covariant and yields complete informa- 
tion about each of the coefficients in the expansion, indicating that all of the terms in a large 
mass DeWitt-Schwinger expansion proportional to inverse, even powers of m have been de- 
termined by the generalized expansion. These terms are constructed out of real curvature 
and electromagnetic field tensors and thus are incapable of yielding the imaginary part of 
G^. Yet the definition of G^ for a complex scalar field, 

GW(x,x') = (0\{<f>(x),<f>*(x')}\0), (64) 

shows that it is, in general, a complex quantity. Any method for calculating G^ that does 
not yield both real and imaginary pieces is incomplete in its determination of G^ . DeWitt 
emphasized this limitation existed for the point-splitting procedure by explicitly stating that 
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the expansion in inverse powers of m 2 was incapable of yielding the imaginary part of G^ 
0. Yet this imaginary part is non- vanishing even in flat space with a constant gauge field. 
Classical field theory calculations of the current due to a complex scalar field under these 
conditions is not identically zero |26|. Schwinger's work incorporated these flat space, 
constant field conditions specifically with the result being the non- vanishing Schwinger pair 
creation rate of Eq.(|5|). 

The only way in which the imaginary part ofG« may be obtained in a DeWitt-Schwinger 
expansion is by obtaining the terms proportional to inverse, odd powers of m. It is possible to 
anticipate the form of the terms in the expansion for G*- 1 '. Moving to units where c — h—1 
but G ^ 1, then G^ and curvature tensors such as R^ v are second order quantities with 
units of {length)' 2 while m and A T have units of {length)' 1 . Using such power counting, 
Davies et al. |]27| constructed the stress energy tensor of a conformally invariant scalar 
field in a conformally invariant spacetime. In the present case of the expansion for G^ 
the same procedure may be used. The term proportional to m _1 must be proportional to 
{length)' 3 in order for G^ to remain second order. Since Eq.(|64"D holds in both flat and 
curved spacetimes, the present method should yield a non-vanishing imaginary part when 
f(r) = h{r) = 1 and when the gauge field A T = constant. Under these conditions, the m _1 
term of Ge should be proportional to {eA T ) 3 and {eA T )r~ 2 . This power counting may be 
expanded by allowing the gauge field to be a function of r, with the result that additional 
terms such as {eA' T )r~ l should appear in the expansion. 

The m~ l term in the DeWitt-Schwinger approximation for GE,ren m the present work is 
given in Euclidean space by 
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All of these terms are proportional to odd powers of A T and its derivatives and, due to 
Eg. (p5|) , are thus all imaginary when rotated back to Lorentzian space. All of these terms 
are also proportional to odd powers of the charge e of the quantum field. Upon imposing 
charge conjugation symmetry on ((J) 2 ) = ^G^ 1 ' PS|| , these terms will not contribute to the 
calculation of the vacuum polarization in the spacetime. 

It is remarkable that the two methods yield exactly the same result for the m~ 2 term 
in their respective DeWitt-Schwinger expansions. As pointed out by Gibbons ||29|| , while 
the DeWitt-Schwinger point-splitting expansion is widely used in renormalization theory, 
there is still much that is unknown about the type of series expansion that it really is. The 
nonzero results of Eq.(|65D highlight the major difference between the present method and the 
generalized DeWitt-Schwinger point-splitting expansion. A major limitation of the DeWitt- 
Schwinger point-splitting expansion is that it is an asymptotic expansion limited to inverse 
powers of m 2 . This dictates that it is restricted to the calculation of only the real part of 
which, as stated above, is generally complex. Yet the DeWitt-Schwinger expansion has the 
distinct advantage of being a completely general expansion with its pieces constructed from 
curvature and electromagnetic tensors. Thus, once the point-splitting calculations have been 
performed and renormalization achieved by discarding terms proportional to nonnegative 
powers of m, the remainder is a completely general expression which may then be used to 
study any spacetime for which a metric exists. 

The main limitation of the present method is that it is restricted to spacetimes for which 
the basis functions of the quantum field can be put into the WKB form like that of Eq.(j3§|). 
The wave equation must therefore be separable. Since this has been achieved for the Kerr 
metric [[30[| , the present work could conceivably be extended to Kerr spacetimes. 



The major advantage of the present method over the DeWitt-Schwinger point-splitting 
expansion is that it can yield the imaginary part of G^\ The calculation of the imaginary 
part of G^ is important due to the definition of the current in Eq. ([T3|) , which is non- 
vanishing even in flat space with a constant gauge field. The factor of i in front of the 
expression, along with the subtraction of the complex conjugate of the derivatives of G^ l \ 

24 



show that only the imaginary part of will contribute to the current. The factor of 
e in front of the expression for requires the presence of the odd powers of e in the 
imaginary part of if the current is to be constructed from G^\ The final expression for 
the current must be proportional to even powers of e if the current is not to vanish under 
charge conjugation symmetry f28fl . Work is in progress to calculate the DeWitt-Schwinger 



approximation for the current in two ways; first, directly from Eqs.fl65|) and (0), and second, 
starting from Eq.([0J) and proceeding through the steps outlined in this paper to calculate 
the DeWitt-Schwinger approximations for quantities such as G^ 1 ^ that are required to 
construct 

Much of the previous research in quantum field theory in curved spacetimes dealt with 
the production of virtual particles due to vacuum polarization. Whenever real particle 
production due to an electromagnetic field was to be studied, investigators most often used 
the pair creation rate derived by Schwinger for flat spacetime ||31||. While Schwinger's work 



assumed a constant electric field E, the authors in Ref. []3l| replaced that field with the 



spherically symmetric radial field of the Reissner-Nordstrom spacetime. With the curved 
space results in Eq.(^5]), it should now be possible to determine how the effects of curved 
space will cause the rate of charged particle production to differ from that of the flat space 
Schwinger rate. 

ACKNOWLEDGMENTS 

The author would like to thank William A. Hiscock and Paul R. Anderson for helpful 
discussions concerning this work. This work was supported by a Radford University Seed 
Grant. 



25 



REFERENCES 

* electronic mail address: rherman@runet.edu 

[1] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (Freeman, San Francisco, 
1973). 

[2] N. D. Birrell and P. C. W. Davies, Quantum fields in curved space (Cambridge University 
Press, Cambridge, England, 1982). 

[3] S. S. Schweber, An Introduction to Relativistic Quantum Field Theory (Row, Peter- 
son, and Company, New York, 1961). See chapter 16 for a discussion of the degree of 
divergences in QED. 

[4] J. D. Bjorken and S. D. Drell, Relativistic Quantum Mechanics (McGraw-Hill, New 
York, 1964). 

[5] J. D. Bjorken and S. D. Drell, Relativistic Quantum Fields (McGraw-Hill, New York, 
1965). 

[6] B. S. DeWitt, Phys. Rep. 19C(1), 297 (1975). 
[7] J. S. Schwinger, Phys. Rev. 82 664 (1951). 

[8] J. S. Schwinger, Phys. Rev. 82 664 (1951). Section VI contains the relevant perturbation 
theory. 

[9] B. S. DeWitt, The Dynamical Theory of Groups and Fields (Gordon and Breach, New 
York, 1965). 

[10] S. M. Christensen, Phys. Rev. D 14, 2490 (1976). 

[11] S. M. Christensen, Ph. D. thesis, University of Texas at Austin, 1975. 

[12] P. Candelas, Phys. Rev. D 21, 2185 (1980). 

[13] V. P. Frolov, Phys. Rev. D 26, 954 (1982). 

26 



[14] V. P. Frolov, in Quantum Gravity, proceedings of the Second Seminar on Quantum 
Gravity, Moscow, USSR, edited by M. Markov and P. C. West (Plenum, New York, 

1984) . 

[15] P. Candelas and K. W. Howard, Phys. Rev. D 21, 1618 (1984). 
[16] P. Candelas and B. P. Jensen, Phys. Rev. D 33, 1596 (1986). 
[17] P. R. Anderson, Phys. Rev. D 41, 1152 (1990). 

[18] K. W. Howard and P. Candelas, Phys. Rev. Lett. 53, 403 (1984); K. W. Howard, Phys. 
Rev. D 30, 2532 (1984). 

[19] V. P. Frolov and A. I. Zel'nikov, Phys. Rev. D 29, 1057 (1984). 

[20] M. R. Brown, A. C. Ottewill, and D. N. Page, Phys. Rev. D 33, 2840 (1986). 

[21] P. R. Anderson, W. A. Hiscock, and D. A. Samuel, Phys. Rev. D 51, 4337 (1995). 

[22] Rhett Herman, Ph. D. thesis, Montana State University, 1996. 

[23] R. Herman and W. A. Hiscock, unpublished. 

[24] Rhett Herman and William A. Hiscock, Phys. Rev. D 53, 3285 (1996). 

[25] E. T. Whittaker and G. N. Watson, A Course of Modern Analysis (Cambridge University 
Press, Cambridge, England, 1927). 

[26] L. H. Ryder, Quantum Field Theory (Cambridge University Press, Cambridge, England, 

1985) . 

[27] P. C. W. Davies, S. A. Fulling, S. M. Christensen, T. S. Bunch, Ann. Phys. 109, 108 
(1977). 

[28] Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, and E. Mottola, Phys. Rev. D 45, 
4659 (1992). 

27 



[29] G. W. Gibbons, in General Relativity: An Einstein Centenary Survey, edited by S. W. 
Hawking and W. Israel (Cambridge University Press, Cambridge, England, 1979). 

[30] B. Carter, Phys. Rev. 174, 1559 (1968). 

[31] See, for example, I. D. Novikov and A. A. Starobinskii, Sov. Phys. JETP 51, 1 (1980), 
and W. A. Hiscock and L. Weems, Phys. Rev. D 41, 1142 (1990). 



28 



